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The filled skutterudite superconductor PrPt4Gei2 was studied in muon-spin rotation (/^SR), spe- 
cific heat, and electrical resistivity experiments. The continuous increase of the superfluid density 
with decreasing temperature and the dependence of the magnetic penetration depth A on the mag- 
netic field obtained by means of £iSR, as well as the observation of a T 3 dependence of the electronic 
specific heat establish the presence of point-like nodes in the superconducting energy gap. The 
energy gap was found to be well described by A = Ao|fca: ±ik y \ or A = Ao(l — ky) functional forms, 
similar to that obtained for another skutterudite superconductor, PrOs4Sbi2. The gap to T c ratios 
were estimated to be Ao/fcsT c = 2.68(5) and 2.29(5), respectively. 

PACS numbers: 76.75.-K, 74.70.Dd, 74.25.Ha 
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The filled skutterudite compounds RM^Pnyi [R = A 
rare-earth or alkaline-earth, M = Fe, Ru, Os, and Pn = 
P, As, Sb) have attracted much attention in the recent 
years. Depending on the composition they may undergo 
metal-insulator transitions show conventional or un- 
conventional superconductivity [3, 0, Hj], become mag- 
netic SB, @|, semiconducting or exhibit the Kondo 
effect [9[ etc. One of the most interesting discoveries 
was superconducting Pr-fillcd skutterudites, since within 
the traditional mechanism of s— wave Cooper pairing the 
Pr magnetism would destroy superconductivity. The 
first compound demonstrating such unusual behavior was 
PrOs4Sbi2 0. The Cooper-pairing mechanism and the 
corresponding symmetry of the superconducting energy 
gap in PrOs4Sbi2 are still under debate. Studies of the 
thermal conductivity suggested the presence of two dis- 
tinct superconducting phases [1(3] • Penetration depth ex- 
periments were pointing either to a possible nodal struc- 
ture of the gap [ll[ or to a gap without nodes A 
number of experimental techniques such as scanning tun- 
neling microscopy thermal conductivity and 
specific heat measurements (l4j are in agreement with a 
fully developed isotropic s— wave superconducting gap. 

Recently, a new Pr-containing skutterudite supercon- 
ductor PrPt4Gei2 was discovered [15j. Specific heat ex- 
periments reveal strongly coupled superconductivity with 
a transition temperature T c ~ 7.9 K, factor of 4 higher 
than that in PrOs4Sbi2 Q, thus making PrPt4Gei2 more 
accessible for e.g. spectroscopic studies. Here, we report 
on a study of PrPt4Gci2 by means of muon-spin rota- 
tion (/zSR), specific heat, and electrical resistivity. The 
linear increase of the superfluid density (p s ) with decreas- 
ing temperature, its dependence on the magnetic field, as 
well as the observation of a T 3 dependence of the elec- 
tronic specific heat document presence of the point-like 
nodes in the energy gap. The temperature dependence 
of p s was analyzed with various gap models suggested 
previously for PrOs4Sbi2 and found to be well described 
by two models of axial symmetry with point-like nodes: 



Ao|fca; i ik y \ or A 



Ao(l - K). 



The sample preparation procedure of PrPt4Gei2 is de- 
scribed in [l5j. The transverse field (TF) ^iSR experi- 
ments were performed at the 7rM3 beam line at the Paul 
Scherrer Institute (Villigen Switzerland). The sample 
was field-cooled from above T c down to 1.5 K and mea- 
sured as a function of temperature in a series of fields 
ranging from 35 mT to 640 mT. Additional experiments 
down to T ~ 0.03 K were performed at 75 mT. Typical 
counting statistics were ~ 7 x 10 6 positron events per 
each particular data point. Electrical resistivity R(T, H) 
and specific heat C(T) down to 0.4 K were measured in a 
commercial system (PPMS, Quantum Design) using an 
AC bridge (LR-700, Linear Research) and the HC option 
of the PPMS, respectively. 

Figure [1] shows the electronic specific heat C e \/T vs 
T 2 . The phonon contribution oc T 3 was subtracted us- 
ing a Debye temperature 0d = 189 K. The normal state 
electronic term is 7^ = 63(2) mJmol _1 K -2 . At the 
lowest temperatures an upturn from a nuclear Schottky 
contribution of Pr becomes visible (C nuc i oc T~ 2 ). Af- 
ter its subtraction we find C e i(T) = j'T + ijT 3 with 
7' « 1.33 mJ mol -1 K~ 2 originating from a minor metal- 
lic impurity phase. For T/T c < 0.2 the clear T 3 de- 
pendence of the superconducting state electronic term 
(r) = 2.93(2) mJmol -1 K~ 4 ) suggests that the gap func- 
tion of PrPt4Gei2 has point-like nodes [T3|. The further 
detailed investigation of the energy gap and the order 
parameter symmetry was performed in TF /iSR experi- 
ments. 

The TF /iSR data were analyzed as follows. The spa- 
tial magnetic field distribution within the flux line lattice 
(FLL) was calculated in a standard way via: 



B(r) = <B)]Texp(- 

G 



iG ■ r)Ba(X, £, b). 



(1) 



Here, (B) is the average magnetic field inside the super- 
conductor, £ the coherence length, b = (B) /B c2 the re- 
duced field, r the vector coordinate in a plane perpendic- 
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FIG. 1: (Color online) Inset: specific heat C — C p h of 
PrPt4Gei2 corrected for the phonon term. Main panel: data 
below 1.4 K and data corrected for the nuclear Schottky con- 
tribution of Pr (see text). 
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FIG. 2: (Color online) The magnetic field distribution P(B) 
at T = 1.55 K and applied fields B a = 75, 200, and 400 mT 
obtained by means of fast Fourier transform. The solid lines 
are fits to the data (see text for details). 



ular to the applied field, G = 47r/v / 3a(m\/3/2, n + to/2) 
the reciprocal lattice vector of the hexagonal FLL, a the 
intervortex distance, and m and n are integer numbers. 
The Fourier components Bq were obtained by minimiz- 
ing the Ginzburg-Landau (GL) free energy using the nu- 
merical algorithm described in Ref. Il7|. The ^tSR time- 
spectra were fit to a theoretical polarization function P(t) 
by assuming the internal field distribution Pid(B) ob- 
tained from Eq. |T]) and accounting for the FLL disorder 
by multiplying P;d(i?) to a Gaussian function [Hi ]: 

P(t)=Ae i * /V (,T ' +<T — )t2 l 2 P iA {B)e l ^ m dB. (2) 



Here A and <f> are the initial asymmetry and the phase of 
the muon spin ensemble, a g is a parameter related to FLL 
disorder [l9( and o nm the nuclear moment contribution 
measured at T > T c . For a detailed description of the 
fitting procedure we refer to [2fj| . 

Figure [5] shows representative internal field distri- 
butions P{B) obtained from the measured /xSR time- 
spectra by performing the fast Fourier transform. The 
solid lines correspond to the fits by means of the above de- 
scribed procedure. P(B)'s have asymmetric shape as ex- 
pected for a field distribution within well arranged FLL. 
The small peak in the vicinity of the applied field is due to 
muons stopped outside of the sample (~2-3%). Above T c 
an additional magnetic depolarization is observed, which 
increases with increasing field. The measurements at zero 
field and above T c revealed that this magnetic depolar- 
ization is temperature independent. 

Magnetic penetration depth is defined in the London 
limit of low magnetic fields as a measure of superfluid 
density p s oc 1/A 2 . The A obtained in the vortex state is 
often called an effective magnetic penetration depth A e ff 
pll. I22I] , since it might be field dependent. For a classical 



BCS superconductor with isotropic s-wave gap A ff = A 
[22] . The fit of Eq. © to the single /iSR time-spectra 
does not allow to obtain independently A and £, since 
they strongly correlate [2(| . In order to get rid of this cor- 
relation one needs to perform a simultaneous fit of several 
spectra, measured at the same temperature but differ- 
ent mag netic fields, with A and £ as common parameters 
[Til |2Q| . B C 2 = $0 /27r£ 2 as a function of temperature, 
obtained after such fit, is presented in Fig. [3k,. The fit 
was performed by combining the data points measured at 
the applied fields B a = 35 and 400 mT. Figures [3)3 and c 
show the corresponding correlation plots [x 2 (A,£)]. The 
normalized \ 2 was obtained for A and £ varied within 
10% around their optimal fit values [Ao = 114(4) nm 
and £0 = 18.1(8) nm]. Strong correlations between A 
and £ as well as a field dependent slope of the correla- 
tion curves are obvious. Figure^ implies, however, that 
B c i (T) obtained in the simultaneous fit is approximately 
40% smaller than that measured directly in specific heat 
[TH ] and resistivity experiments. This suggests that the 
assumption of field-independent A may not be valid for 
PrPt4Gei2 and A depends on field. Field dependence of 
A is expected for a superconductor with nodes in the gap. 
A field induces excitations at the gap nodes due to nonlo- 
cal and nonlinear effects, thus reducing superconducting 



carrier concentration n s oc 1/A 2 [2l|. We believe how- 
ever that GL theory developed for superconductors with 
isotropic s-wave gap is approximately valid in the present 
case. In order to avoid the correlation between A and £ 
we fixed the values of £(T) = v /$o/27rS c2 (T) by using 
the B C 2(T) curve obtained in resistivity measurements 
(see Fig.lSt) HH- 

The GL theory allows to extend the London definition 
of superfluid density for a finite applied field [23| : 
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FIG. 3: (Color online) (a) Temperature dependence of the up- 
per critical field B C 2 obtained in specific heat (squares) and 
resistivity (dots) experiments compared with results of the 
simultaneous fit of pSK time spectra at 35 and 400 mT (tri- 
angles). Panels (b) and (c) show correlation plots X 2 (-V£) f° r 
B a = 400 mT and 50 mT, respectively. 




FIG. 4: (Color online) Field dependence of the superfluid 
density p s oc (|t/>| 2 }/A 2 . The solid line is the power law fit 
p a (T)/p a (T = 0) = 1 - (B/B c2 ) a with a ~ 0.25. The dashed 
line represents the case when A does not depend on the mag- 
netic field. The inset shows A(B/B c2 ) at T = 1.55 K. 



Ps cx(|^r)| 2 )/A 2 c(l-fe)/A 2 



(3) 



Here, tp( r ) is the GL order parameter and (...) means 
averaging over the unit cell of the FLL. 

Equation ([3]) implies that in an isotropic s-wave 
superconductor the mean superfluid density decreases 
with increasing field as (|-0(r)| 2 ) (the dashed line in 
Fig. dj. Obviously, the experimental p s decreases much 
stronger with increasing magnetic field than for a field- 
independent A, thus suggesting the presence of nodes in 
the gap. The power law fit p s (T)/p s (T = 0) = 1 - b a 
results in a = 0.26(2) (solid line in Fig. |4j) . 

Figure shows p s , normalized to its value at T = 
30 mK, as a function of the reduced temperature t = 




FIG. 5: (a) Temperature dependence of p s oc (\ip(r)\ 2 }/\ 2 
normalized to its value at T = 30 mK. The lines represent 
the fits by using 5 different gap models. The inset shows the 
low temperature region between T — and T = 0.3T C . The 
panels from (b) to (f ) represent the gap functions used for the 
analysis of p s (T) data (see Table U for details). 



T/T c . Surprisingly, the t dependence of p s is identical 
for all applied fields. Only the data points measured at 
35 mT deviate slightly. This may be due to larger vortex 
disorder and, consequently, enlarged systematic errors at 
such a low field. The inset of Fig. [5^ shows p s measured 
at 75 mT for T < 0.3T C . It is obvious that p s does not 
saturate but increases continuously with decreasing tem- 
perature as expected for a superconductor with nodes in 
the gap. Note that the presence of nodes is also consis- 
tent with the conclusion drawn from the analysis of p s (b) 
(see Fig. [4] and the discussion above). 

We are not aware of any work reporting the possi- 
ble order parameter symmetry in PrPt4Gei2. This is- 
sue was thoroughly studied for another representative 
of the skutterudite family PrOs4Sbi2 
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Conse- 
quently, the temperature dependence of the superfluid 
density was analyzed by assuming the following gap sym- 
metries: (1) the isotropic gap A(6, 4>) = Ao (Fig.^b), (2) 
A(0,0) = A \k x ±iky\ = A osin# (Fig. [5b), (3) A(0,0) = 
A (l -k£- = A (l - sin 4 (9 cos 4 - sin 4 e»sm 4 ^) 



(Fig. EH), (4) A( 



A (l-fc 4 ) = A (l-sin 4 6»sin 4 



4 



TABLE I: Summary of the gap analysis of p 3 (T) for 
PrPt4Gei2 (see text for details). 
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1.33(4) 1.95(5) 


3.79 
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1.82(4) 2.68(5) 
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1.56(4) 2.29(5) 


1.12 


Ao|fcy 
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3.32(20) 4.88(30) 


8.85 


A (l-fe*-fc*) 
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2.61(10) 3.84(14) 


3.29 



(Fig. Eh), and (5) A(6,<p) = A k y = A | sinflsin^ 
(Fig. Ef)- Here Ao is the maximum value of the gap 
at T = 0, and 9 and <f> are the polar and the az- 
imuthal coordinates in the fc-space. The temperature 
dependence of the gap A(#, </>, i) = A(9,<fi)5(t) was as- 
sumed to follow the weak-coupling BCS prediction 5(t) = 
tanh{1.82[1.018(i- l) 51 ]} El- Note, 5(i) is practically 
independent of gap model [26|. Although the gaps de- 
scribed by (2) and (3) have similar shape (see Figs. Efc 
and d), the A = A (l — k y ) one corresponds to triplet 
pairing [24| which is observed in PrOs4Sbi2 [!]• 

The temperature dependence of the superfluid density 
was calculated within the local (London) approximation 
by using the equations [26| : 



p aa = 1- 



Pc 



AttT 



2itT 



sin 



cos- 



cosh 



2T 



cos 2 9 cos 2 6 cosh 



ded(j)d9 
(4) 

ded<pd9 
(5) 

Here, pa (i = a, b, or c) is the component of the superfluid 
density along ith principal axis. Since our experiments 
were performed on a polycrystallinc sample, we used the 
powder average of the superfluid density: 



2T 



Ps = (^PaaPbb + y/PaaPcc + y/PccPbb)/^- 



(6) 



The fit was made for the only free parameter Ao. The 
results of the fit are presented in Fig. Eh and Table [T] 
The goodness of fit was checked by using the \ 2 cr i _ 
tcrium. The values of \ 2 normalized to that obtained 
for the best Aoj/Ca; ± ik y \ model are given in the last col- 
umn of Table U ft is obvious that the fits using gaps 
of the form: A = A (l - k% - k*), A = A k y , and 
A = Ao are very poor. The smallest \ 2 was obtained for 
the A = Ao|fca; ± ik y \ model, while the A = Ao(l — k y ) 
model results in ~12% larger x 2 value. The gap-to- T c ra- 
tios were estimated to be Aq/ ksT c = 2.68(5) and 2.29(5) 
for the Ao|fez ± ik y \ and Ao(l — k y ) gap functions, re- 
spectively. Both ratios are higher than the BCS weak 
coupling value 1.76, thus suggesting that PrPt4Gci2 is 



a moderately strong coupled superconductor. They are 
also relativel y c lose to A/ksT c = 2.35 obtained from 
specific heat [151 ] . 

To conclude, the PrPt4Gei2 skutterudite superconduc- 
tor was studied in pSR, specific heat, and electric resistiv- 
ity experiments. For T/T c < 0.2 a clear T 3 dependence 
of the electronic specific heat term gives evidence that 
the gap function of PrPt4Gei2 has point nodes. Analysis 
of the /xSR data was performed by exact minimization of 
the Ginzburg-Landau free energy. The dependence of the 
superfluid density (p s ) on the magnetic field was found 
to follow p s cx 1 — (B j B C 2) 25 . The temperature depen- 
dence of p s (T) was found to be well described by two 
models with axial symmetry with point-like nodes in the 
gap: A = Aol^ ± ik y \ and A = A (l — k y ) despite of 
cubic symmetry of crystal structure. The maximum gap- 
to- T c ratios were estimated to be A /kBT c = 2.68(5) and 
2.29(5), respectively, in agreement with 2.35 reported in 
Ref. [l5j based on specific heat results. 

This work was performed at the Swiss Muon Source 
(S/xS), Paul Scherrer Institut (PSI, Switzerland). We 
would like to acknowledge E. H. Brandt for valuable dis- 
cussions. 
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